Sierpinski Gaskets for Logic Functions Representation 
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Abstract 

This paper introduces a new approach to represent logic 
functions in the form of Sierpinski Gaskets. The structure 
of the gasket allows to manipulate with the corresponding 
logic expression using recursive essence of fractals. Thus, 
the Sierpinski gasket's pattern has myriad useful proper- 
ties which can enhance practical features of other graphic 
representations like decision diagrams. We have covered 
possible applications of Sierpinski gaskets in logic design 
and justified our assumptions in logic function minimiza- 
tion (both Boolean and multiple-valued cases). The exper- 
imental results on benchmarks with advances in the novel 
structure are considered as well. 



1 Introduction 

Over the years, many important problems in digital cir- 
cuits synthesis have been approached using graph-based 
data structures: decision trees, decision diagrams . . . etc. 
The superior example is Binary Decision Diagram (BDD) 
which has become the advanced structure in VLSI CAD 
for representation and manipulations of logic functions [|lj. 
The applications of BDD techniques cover practically all 
stages of digital circuits design from the initial represen- 
tation of Boolean functions [0] to the synthesis of the de- 
sirable circuit [H], But still originally introduced BDDs 
cannot be directly applied to synthesize circuits derived 
from AND/EXOR expressions. Such circuits are demanded 
by VLSI CAD, because of economical implementation (in 
terms of gates and interconnections) and high testability. 
This is particularly efficient for error control and arithmetic 
circuits, encrypting and coding schemes. In this paper, 
we address the problem of representation of the general 
class of AND/EXOR expressions, named Exclusive Sum- 
Of-Products (ESOP). Following the above problem, we in- 
troduce another graph-based structure popular in fractal 
theory to represent ESOP expressions and their multiple- 
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valued counterparts. 

One is familiar with the useful properties of Pascal tri- 
angle which represents binomial coefficients. Another ad- 
vanced mathematical structure, which is used frequently in 
fractal theory, is Sierpinski gasket [§], It should be men- 
tioned that Sierpinski gasket is Pascal triangle modulo two - 
the EXOR operation is used instead of addition [j^|. Sierpin- 
ski gasket also referred to as Sierpinski triangle has a num- 
ber of interesting properties that W.Sierpinski discussed in 
his paper [pp. Areas of application of Sierpinski fractals 
are graphic design (Figure [l] shows (a) Sierpinski gasket 
and (b) Sierpinski pyramid as the samples of fractals imple- 
mentation) [n], telecommunication (the antenna of cellular 



phone has fractal square structure), tautology mapping [ 1 1 1, 
symmetry handling in logic design [Q, |l2|] . . . etc. In this 
paper, we propose a new approach to utilize the structure of 
Sierpinski gasket to represent ESOP expressions. Prelimi- 
nary case study results have been published in [Bp. We also 
generalize this approach to multiple-valued functions. 

The rest of the paper is organized as follows. In Sec- 
tion ^, we collect necessary definitions and provide basic 
terminology. Section ^] introduces the structure of Sierpin- 
ski gasket and describes the techniques of constructing Sier- 
pinski gasket using different function descriptions. The cor- 
respondence between Sierpinski gaskets and logic expres- 
sions is outlined in this Section. Section [| is dedicated 
to a minimization algorithm and experimental results on 
Boolean and multiple-valued functions. Finally, Section [| 
concludes the paper. 

2 Basic Concepts 



We consider a multiple-valued function f: M" — » M 
over the variable set X = {xi, ■ ■ ■ , x n }, where 
M={0, 1, . . . , r — 1}, n is the number of r-valued vari- 
ables. f\ Xl=v = f(xi, . . .,Xi-x,v,x i+ i, . . . ,x n ) is called 
a cofactor of /, when Xi takes value v 6 {0, 1, . . . i — 1}. 

In the Galois field representation of MVL functions, p 
complements of a p— valued variable are usually consid- 




Figure 1. Sierpinski gasket and pyramid 

ered. They are defined by l ~x = x + i, i = 1, ... ,p — 1, 
where '+' denotes the addition in the considered Galois 
field. Therefore, for functions over GF(4), we consider four 
literals 30 ^ 30 j 30 ^ 3C for each variable x. 

For simplicity, the presentation in this paper is restricted 
to Boolean and 4- valued functions. However, our method 
can be applied to logic functions of any radix. Let us review 
the representation forms of logic functions. 

2.1 Boolean functions 

It has been shown in [^J that restricted classes of 
AND/EXOR expressions, like Fixed Polarity Reed-Muller 
(FPRM) expressions, can be represented efficiently. 

Definition 1. FPRM is an exclusive OR of AND product 
terms, where each variable appears either complemented 
or uncomplemented, but not both. 

Based on Sierpinski gasket we are able to represent gen- 
eral class of AND/EXOR expressions - ESOR 

Definition 2. ESOP is EXOR sum of product terms with 
complemented and uncomplemented variables. 

Observe that in ESOP expression a variable may appear 
in a product term in complemented or uncomplemented 
forms. In case of ESOP, the polarity vector [cqC\ . . . c„] 
is associated with given function /, where the entry = 
{0, 1 , d} is if Xi appears always in uncomplemented form, 
the entry is 1 for the variable in complemented form and is 
d if Xi is used in complemented as well as uncomplemented 
forms (mixed mode) [Q] . Note that different polarity vectors 
for the same function cause different cost characteristics of 
AND/EXOR expressions (number of product terms, num- 
ber of literals) [§]. 

Example 1. For the function f = x\ V X3 ■ X2 with the 
truth vector [10001111], FPRM with polarity vector [111] 
is f = 1 © X\ ffi X1X2X3 ( 3 product terms), and ESOP with 



polarity vector [dll] is f = x\ (BX1X2X3 (2 product terms). 
The effect in number of product terms has been obtained for 
different polarity vectors. 

2.2 Multiple-valued functions 

The fixed polarity expansion for 4-valued function can 
be introduced. 

Definition 3. (D. Green /Q/j Reed-Muller expansion of 4- 
valued function takes the form 

/ = Co + cix + c 2 x 2 + c 3 x 3 , 

where operations are fulfilled in GF(4). 

In our approach, we extend the representation of 4- 
valued functions using mixed polarities: 

/ = c + ci-~x + c 2 1 ~x 2 + c 3 1 'x 3 , 

/, 2 — 1 2 — 2 1 2 — 3 
= c + ci x + c 2 x + c 3 x , 

/ = c + ci^x + c 2 3 ~x 2 + c 3 3 ~x 3 . 

All possible combinations of four expansions specified 
above can be represented using Sierpinski pyramid. 

3 Sierpinski Gasket 

Definition 4. Sierpinski gasket in our study is a recursively 
connected graph with the vertex set and the edge set, where: 

(i) Each vertex is labeled by a product term p assigned, as 

a decision term. Each vertex associated with product 
term p has at least two connected vertexes p\ and p 2 
correspondingly which are related using the following 
triangle rule: 

P®Pl®P2=0. (1) 

(ii) Each triangle contains three sub-triangles except the 
primitive one. The biggest equilateral triangle in Sier- 
pinski gasket encloses three equilateral triangles with 
the base one third of the size of the original, and with 
hole in the center. Replacement of each of those sub- 
triangles gives three more sub-triangles, to obtain the 
gasket as nine triangles, each with a base of length 1/9. 
The repeating of this replacement n — 1 times produces 
the final gasket. 

The vertexes of Sierpinski gasket ar e ad dressed accord- 
ing to a coordinate convention (Section 3.1 ). 



Example 2. The gasket for the function f from Example [7] 
is given in Figure || 
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Figure 2. Sierpinski gasket for Example Q 




Figure 3. Sierpinski gasket and the coeffi- 
cients of positive polarity and negative po- 
larity Reed-Muller expressions 



As well as other fractals, Sierpinski gasket can be defined 
as a recursive procedure. We propose to use the following 
recursive procedure to build Sierpinski Gasket (SG) for the 
function / (see Definition ||for details): 

SG(f) = SG(f {x=0 ) U SG(f ]x=1 ) U SG(f ]x=d ) 

Another technique for building Sierpinski gasket is to 
use the truth vector of the function /: 

Step 1. Start building by arranging elements of the truth 
vector at the base level {level = 0). 

Step 2. Select all pairs of elements, starting from the left 
side, EXOR them, and than write them down at the 
upper level. 

Step 3. Repeat the second step until the level contains one 
element only. 



binary code of the product term of variables without 
negation, and the second j is equivalent to the binary 
code of the product term with complemented variables. 

Lemma 1. The sum of both coordinates i and j is 
equal to level number of the Sierpinski gasket. 

Theorem 1. The total memory size (in elements) that 
is needed to represent a function f using Sierpinski 
gasket and the coordinate axes described above is 
equal to ^ 4p 

Theorem 2. The number ofvertexes in Sierpinski gas- 
ket is equal to 3™ for the n-variable function f. 

Corollary 1. The number of empty elements in Sier- 
pinski gasket is equal to 4; ^ 3". 



Example 3. Consider a two variables function f specified 
by truth vector [abed]. Figure ^a) shows the base level 
(truth vector) for constructing Sierpinski gasket. The coef- 
ficients {a © b © c © d, b © d, c © d, d} of positive polarity 
Reed-Muller expression for given function are represented 
by the right side of the triangle (Figure $b)), coefficients 
{a © b © c © d, a © c, a © b, a} of negative polarity Reed- 
Muller expression are located on the left side of the triangle 
( Figure %c)). 



3.1 Properties 

To manipulate with Sierpinski gasket as regular fractal 
structure, we need to apply a coordinate convention. Two 
options can be analyzed. 

Coordinate convention 1. Each element (vertex) of the 
Sierpinski gasket can be addressed as two-dimentional 
array where the first i index corresponds to the 



Coordinate convention 2. Another coordinate convention 
assumes that each element of Sierpinski gasket has n- 
dimensional coordinates. For Boolean functions, we 
have to encode three states of the variable x: d (the 
variable is not present in an expression), x (variable 
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Figure 4. Different coordinate systems 



with complement) and x (variable without comple- 
ment). The following values are assigned for variable's 
states: '1' ford, '2' for 3; and '3' forx. Such encoding 
has been chosen to provide the similarity to the triangle 
rule (Equation |lj): 1 © 2 3 = 0. 

Theorem 3. The total memory size (in elements) 
needed to represent a n-variable Boolean function f 
using Sierpinski gasket and the n-dimensional coordi- 
nates is 3 n . 

These coordinate systems are shown in Figure [|. To min- 
imize the space complexity, we have selected the second 
coordinate convention to manipulate with Sierpinski gas- 
ket (Boolean function) and Sierpinski pyramid (4-valued 
function). Recall that coordinate convention assumes that 
each element of Sierpinski gasket can be addressed using 
ri-dimensional coordinates. 

3.2 Applications 

3.2.1 Logic Expansions 

Generally, a decomposition of a Boolean function / with 
respect to arbitrary variable x, or, in other words, an expan- 
sion of / given x, can be represented by the formula: 

f = (C -x-f )°(C 1 -x-f 1 )o(C d -f d ), (2) 

where C — [CoCiCd] is a vector of coefficients, and symbol 
< o > denotes a logical operation, EXOR for AND/EXOR 
expressions. Below we will utilize the following expan- 
sions derived from possible combinations of coefficients 
Co,Ci,Cd- Shannon (S) expansion (C = [110]); positive 
Davio (pD) expansion (C = [Oil]), and negative Davio 
(nD) expansion (C = [101]). The expansion rules and their 
multiple-valued counterparts in GF(4) are given in Table |} 
The main difference between the proposed structure of 
Sierpinski gasket and the original one (Figure |(a)) is that 
the novel structure reflects the properties of the whole set 
of expansions. The simple manipulations with Sierpinski 
gasket yield to the rules with triangles (Figure |[). 

3.2.2 ESOP 

Sierpinski gasket contains complete information to recon- 
struct the function in the form of AND/EXOR expression. 
Thus, all coefficients of ESOP expression can be obtained 
from Sierpinski gasket. And vice versa, ESOP expression 
can be uniquely represented by Sierpinski gasket. 

Each vertex of Sierpinski gasket has the equivalent in 
the form of logic expression (product term): (i) the ones as- 
signed to vertexes of Sierpinski gasket correspond to prod- 
uct terms of ESOP expression; (ii) zeros and d values are 
considered for manipulations only. Final expressions can 



Table 1. Shannon and Davio expansions and 
their analogues in GF(4) 



Type 


Rule of Expansion 


S 


/ 


— x ' f\x=0 © x ' f\x=l 


pD 


/ 


= f\x=0 © x ■ (/|x=0 ffi f\x=l) 


nD 


/ 




4-S 


f 
J 


— J0\X) J\ x= q -f- J1(X) J\ x —i-t- 
M x ) ■ f\x=2 + M x ) ■ f\x=3> 


4-pD 


/ 


= f\x=0 + x ■ (f\x=l + 3f\x=2 + 2/| x= 3) 
+X 2 -(/|*=l+2/|,=2+3/|^3) 

+Z 3 ' U\x=0 + f\x=l + f\x=2 + f\x=s) 


1-4-nD 


/ 


= f\x=l + 1 ~ x ' C/>=0 + 1f\x=2 + 3f\x=z) 
+ 1 ~ x ' 1 ■ {f\ x =0 + 3 f\x=2 + 2 /|iC=3) 
+ 1_ Z 3 ■ (f\x=Q + f\x=l + f\x=2 + /|a=3) 


2-4-nD 


/ 


= f\x=2 + 2 ~ x ■ (3/| a =0 + 2/|as=l + f\x=3) 
+ 2 ~X 2 ■ (2/| x=0 + 3/| a=1 + /| a= 3) 

+ 2 ~ x3 ■ (f\ m =0 + f\x=l + f\x=2 + f\ x =a) 


3-4-nD 


/ 


= f\x=3 + 3 ~ x ■ (2/| a =0 + 3/|j. = i + f\ x=2 ) 

+ z -x 2 ■ (3/ N=0 + 2/ k=1 + f ]x=2 ) 

+ 3 ~X 3 ■ {J\x=0 + f\x=l + f\x=2 + f\x=3) 



be formed using EXOR operation applied for product terms. 
We interpret d as an unspecified coefficient of ESOP expres- 
sion. For arbitrary function / we can build a set of Sierpin- 
ski gaskets with different costs (number of product terms, 
number of literals). Value d allows us to manipulate with 
the set of gaskets, for example, to find the gasket and re- 
ferred ESOP expression with the minimal cost. The manip- 
ulation rules which can lead to extension and simplification 
of the logic expression are presented in Figure ^|(a). 

Example 4. The gasket from Figure ^contains the follow- 
ing ones coordinates: (1,2,2), (3,1,3) and (3,3,2). These 
coordinates produce ESOP expression: W2 ■ X3 © x\ ■ 23 © 

Xi- X 2 - ^3- 

3.2.3 Multiple-valued Functions 

The approach that is presented above has been generalized 
for multiple-valued functions. Figure 0(b) gives a graphi- 
cal representation of Sierpinski gasket for 4-valued function 
(fixed polarity 4 — pD in GF(4)). One can conclude that 
this graph denotes Sierpinski pyramid. Another pyramid 
of a Sierpinski kind will serve as a representation of mixed 
polarity expression for 4-valued function (Figure 0(c)). Ma- 
nipulation rules are given in Figure ^(b). We have consid- 
ered 4-valued functions, hence this approach can be gener- 
alized for functions with any radix. 
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Figure 5. Representation of expansions using 
Sierpinski gasket 



4 Minimization Algorithm and Experiments 

The sketch of the Sierpinski — GFSOP algorithm is 
depicted below. 



4.1 Algorithm Sierpinski 
{Sierpinski — GFSOP) 



ESOP 



1. Generate non-zero vertexes of Sierpinski 
Gasket (Sierpinski Pyramid); 

2 . Apply simplification or extension 
rules to existing Sierpinski Gasket 
(Sierpinski Pyramid) ; 

3. Calculate the cost function (number of 
product terms or literals); 

4 . Choose backtracking if no improvements . 

Note that this algorithm manipulates with non-zero co- 
efficients only applying rules of fractal geometry. 

Example 5. The steps of ESOP minimization for the func- 
tion f = xx V afi • 3?2 with truth vector [1011] are given 



(a) I II 

■| Extension d Extension -| 

/ \ Simplification / \ / \ SimDlification / \ 

d d 1 1 1 

(b) 

1 Extension Q 



Simplification , 

d d — 1 



Simplification 

d J d 1 




Figure 6. Manipulation rules for 4-valued func- 
tions represented by Sierpinski gasket 



in Figure ^. Finally, the gasket contains the following ones 
coordinates: (2,3) and (1,1). These coordinates produce 
ESOP expression: 1 © X\X2. 

4.2 Complexity 

Minimization algorithm presented here is similar to 
EXORCISM - MV3 developed in [0. It should be 
mentioned that the algorithm EXORCISM - MV3 al- 
lows to minimize only multiple-valued input, two-valued 
output functions. However, in this paper, we are consid- 
ering multiple-valued input, multiple-valued output func- 
tions. Our manipulation technique is based on fractal no- 
tations and graph simplifications. 

The space required to store the entire Sierpinski gasket 
for Boolean function / and ESOP minimization is 0(3") 
memory locations, 0(4") to represent and manipulate with 
fixed polarity expressions in GF(4), and 0(16") - for mixed 
polarity representations of 4-valued functions. 



4.3 Experimental Results 
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Our Sierpinski — ESOP (Sierpinski — GFSOP for 
multiple-valued functions) program in C++ implements the 
above described algorithm to minimize ESOP expressions 
and multiple- valued functions. All the experiments have 
been carried out on a 800MHz Pentium PC with 128Mb 
of memory. 

4.3.1 ESOP minimization 

In the first series of experiments, we have selected sev- 
eral LGSynth93 benchmarks to build Sierpinski gaskets and 
minimize ESOP expressions. Thus, for X0r5 with 5 in- 
puts, the number of ones vertexes in the gasket is 5, and 
therefore there is 5 product terms in ESOP expression. For 
9sym with 9 inputs, the number of ones vertexes is equal to 
84, which is the number of product terms of ESOP expres- 
sion. Table || contains the details of comparison between 
the proposed approach and two techniques: (i) symbolic 
manipulation EXORCISM-MV3 @ and (ii) DD based 
Est /Greedy [f§. In many cases, our program demon- 
strated superior results. Memory allocation for the selected 
set of benchmarks using Sierpinski gaskets is 1 .4 times ef- 
fective than using BDDs (for instance, 14 bytes for X0r5 
and 292 bytes for 9sym instead of 96 bytes and 400 bytes 
correspondingly in case of BDDs). 

4.3.2 Minimization in GF(4) 

In the second series of experiments, we have tested 
Sierpinski — GFSOP on several 4-valued benchmarks 
(Table |]). The 4-valued benchmarks were generated by 
pairing inputs and outputs of MCNC benchmarks. The ex- 
perimental results demonstrate that Sierpinski — GFSOP 
produces the fewer number of terms and literals against the 
program based on information theoretic measures Info — 
MyPSDRMGF ( see Q f or details), however the program 
Info—MV is extremely faster. 

5 Concluding Remarks 
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We have presented a novel approach for representing 
logic functions using fractals, namely the structures of Sier- 
pinski gasket and Sierpinski pyramid are used. Several re- 
lated issues such as coordinate convention and memory al- 
location have been discussed. Due to the recursive nature 
of Sierpinski gasket, the proposed techniques for construct- 
ing triangles are efficient in both run-time and storage. It 
becomes possible to develop manipulation algorithms, like 
ESOP minimization or minimization of multiple-valued 
functions in GF(4), using representation in the form of Sier- 
pinski gasket. 
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Figure 8. Steps of ESOP minimization using Sierpinski gasket 



Table 2. Experimental results in ESOP minimization 



I/O 


EXORCISM - MV3 llOf 
T/L 


Est/Greedy M 
T/L 


Sierpinski — ESOP 
T/L 


xor5 


5/1 


5/10 




5/5 


rd53 


5/3 


14/57 




17/48 


bw 


5/28 


22/319 


20/- 


26/122 


inc 


7/9 


26/176 


31/- 


32/153 


5xpl 


7/10 


32/170 


47/- 


55/166 


9sym 


9/1 


51/426 




84/449 


apex4 


9/19 


439/6181 




438/3692 


sao2 


10/4 


28/288 


41/- 


55/370 


exlOlO 


10/10 


670/7466 




879/7014 


duke2 


22/29 


78/909 


108/- 


91/793 


x6dn 


39/5 




104/- 


113/1026 


x7dn 


66/15 






594/4610 



Table 3. Comparison of Info-MV [13] and Sierpinski 

Inf^MV PsJ ™ GF Inf(^MV PS1 



GFSOP on 4-valued benchmarks 

KGF Sierpinski - GFSOP 



I/O 


T/L/t 


T/L/t 


T/L/t 


5xpl 


3/5 


165/521/0.04 


142/448/0.39 


140/453/4.53 


9sym 


5/1 


246/800/0.12 


201/776/0.29 


193/764/30.87 


bw 


3/14 


54/144/0.02 


44/132/0.02 


42/140/5.91 


clip 


5/3 


825/3435/4.50 


664/2935/4.51 


606/2802/40.03 


conl 


4/1 


50/138/0.53 


19/50/0.42 


18/48/12.75 


exlOlO 


5/5 


1018/4076/8.06 


997/4985/8.61 


912/3977/145.80 


inc 


4/5 


146/493/0.36 


65/216/0.82 


60/198/15.49 


misexl 


4/4 


48/108/3.6 


15/38/0.02 


15/43/9.76 


rd84 


4/2 


207/656/0.85 


207/656/1.05 


207/656/9.96 


sao2 


5/2 


252/1133/4.8 


96/437/0.24 


81/365/37.06 


squar5 


3/4 


51/135/0.01 


48/128/0.03 


42/112/10.15 



